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Abstract 



This article models the energy cascade singularity process for the incompressible 
Navier-Stokes equations. It is shown that the energy density transport in the 
context of singularity is equivalent to a backwards heat equation. It is then 
shown through a renormalization type argument that this process leads to a 
logical contradiction, which shows that such energy cascade singularity process 
does not exist in the Navier-Stokes equations. This will imply regularity for 
the Navier-Stokes equations as it is shown that the formation of singularity is 
developed through a gradient flux of the energy density always. 



1 Introduction 



The Navier-Stokes equations are supposed to model any kind of fluid flow, in- 
cluding turbulent ones. The essential problem with the equations is the nonlin- 
earity of the convective term. The pressure gradient represents transport and 
the diffusive term represents loss of energy in the system. In other words, to 
understand the equations, one has to understand the nonlinear convection term. 
The Navier-Stokes equations can be written in vector form as follows 

du _, __ . _, _ , . . 

— + u- Vu = vAu - Vp (I) 
at 

where u(x,y, z,t) € M 3 is the velocity field of the fluid and p(x,y,z,t) £ R is 
the scalar pressure field. The diffusive term depends on a constant v > 0. The 
convective term is u ■ Vm. This term can be decomposed according to the vector 
calculus identities 

u -Vu = u) x u + \7(-u ■ u) (2) 
Subsituting this into the Navier-Stokes equations, one has 

du _ , 1 _ . _ . . 

— = — cu x u — \(jp + — u ■ u) + i/Au (3) 

where w — V x u is the vorticity of the fluid flow. This formulation of the 
Navier-Stokes equations is known as the Lamb formulation [T|. In addition we 
assume that the fluid flow is incompressible, that is 

V ■ u = (4) 

The system is assumed to have some smooth and square-integrable initial data 
a° = u(x, y, z, 0) over the whole space R 3 . 

For the well-known facts of the Navier-Stokes system, we refer to [5] 



2 The singularity process 

One of the most celebrated results in the field of Navier-Stokes equations is the 
Beale-Kato-Majda theorem. It states that the if there is a singularity forming 
in the system, then the following holds 

f T 

lim / s)\\ loo ds = oo (5) 

t ^ T Jo 

In other words, if there is a singularity forming, the vorticity is accumulating in 
some region of space quite rapidly. Suppose now that there is indeed a build-up 
of singularity in terms of vorticity, as prescribed in the Beale-Kato-Majda crite- 
rion. As stated earlier, the singularity is forming on the singular set, which can 
be I-dimensional. Nevertheless, as shown in Caffarelli et al.[3], the singular set 



f 



can be covered with a minimal compact set CcR 3 having a smooth surface 
dC. As the Navier-Stokes system is deterministic, the singular set for given 
initial data is definitely some set S C C C M 3 . 

We now therefore assume a situation where the system is forming a singularity 
on some singular set S. The singularity occurs at time T. Let us now take 
ourselves in the time T — e, e > 0, in other words very close to the singularity in 
chronological terms. We cover the singular set S with the set C . As we are very 
close to singularity, the vorticity is already very large on S and its surroundings, 
in particular we choose the minimal covering set C in a way that the enstrophy 
density function is strictly positive in C 

(f> = uj-cu>0 (6) 

which is always possible, as everything is still continuous. In addition, we choose 
the time parameter e small enough such that at T— e the enstrophy density func- 
tion is larger on the singular set than anywhere else on the domain M 3 . 



2.1 The energy cascade 

It is obvious that the velocity field blows up as well when the vorticity blows up, 
see [1], this can be seen from the definition of vorticity: the spatial derivatives 
of the velocity field become infinite by definition at the singularity. This can 
only happen if the velocity field itself becomes infinite at the singularity. As 
stated earlier we consider time T — e with e small enough to ensure that the 
velocity has a global maximum inside the covering set C. 

As we move closer and closer to the singularity moment in terms of time, the 
process of forming of a singularity requires that energy density / = \u ■ u is 
transported from other parts of the system into the singular set, as otherwise 
the the divergence of velocity at the singularity would be impossible. It should 
be noted though that the velocity density is finite inside the set C until the 
moment the singularity occurs. Therefore if we have some minimal compact set 
C in the neighbourhood of S that is embedding the singular set, the net flux 
of energy density through the surface of C must be inwards and thus the flux 
integral must be negative. This means that 

N(t) = I F ■ C < (7) 
Jac 

where F represents the flux vector of energy density. Because the net flux satis- 
fies — oo < N(t) < for all t > T — e, we can always represent the flux of energy 
density through the gradient field of the energy density in the following way. 

By construction, we know that on the surface of the set C the gradient of energy 
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density points to the forming singularity, as we constructed the system in such 
a way that the global energy density maximum is on the singular set at time 
t = T — e. Because we know that -co < N(t) < for all times T - e <t <T 
we can always represent the net flux N(t) with the choice of 



N(t) = & (x(t)V/ • C < (8) 

JdC 

where a(t) > is a scaling function that ensures that the net flux can always be 
represented through the gradient field of the energy density /. For the purpose 
of this study it is however sufficient to freeze the system at t = T — e and study 
the system at that point in time. This means that instead of o~(t) we take just 
some a > so that the net flux at that time instant can be represented through 
the gradient flux. 

Let us now consider what the system does at such time close to singularity. 
The change in energy content of the covering set is given by 



dK(t) 
dt 



This is so, because we have a fixed control volume C in space. So that the 
Reynolds transport theorem is reduced to the equation above. By definition, 
the time derivative of the energy content must be other hand related to the flux 
of energy density 

dK(t) 



dt 



o-Wf-C (10) 

dC 



This tells that the only way that the integral of the energy density within the 
set C can increase is via the flux through the surface, inside the set. As the 
total kinetic energy is decreasing in the system, the increase of kinetic energy in 
some set must brought upon by transfer from other regions of the space. This 
is precisely what the previous equation states. Let us now use the divergence 
theorem. 

/ V • a{t)VfdV = <f> <r(i)V/ • C (11) 

Jc JdC 

As there is a global maximum for / inside the set C, we must have that 

Af < (12) 
on the set C. Combining equations we have 

As the set C is arbitrary, we must have 

% = -<r(t)Af (14) 
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So the energy cascade dynamics is also equivalent to backwards heat equation, 
which is known to behave badly. It can even blow up in finite time, given 
(im)proper conditions. Note that because we are close to singularity in chrono- 
logical terms, the energy density has a global maximum inside the set C . 

3 The Navier-Stokes equations and the singu- 
larity process 

Let us now recall the Lamb formulation of Navier-Stokes equations 
du 1 

— = -L- V(p + -u-u) + vAu (15) 

Where the Lamb vector is defined as 

L = uj x u (16) 

If one now takes divergence, one has the divergence for the Lamb vector 

V-L=-A{p + f) (17) 

So that we can reformulate the backwards heat equation for the energy cascade 
as 

|£=a(i)(V-L + Ap) (18) 
The Lamb vector divergence, or hydrodynamic charge, can be written also as 

V'I = fi'VxiiJ-W'W (19) 

Or in the form 

V • L = -u ■ Au- w ■ uj (20) 
see, [5], So that the energy cascade equation becomes 
df 

-± = cr(t)(Ap-u- Au-uj-uj) (21) 

Let us now take the inner product of velocity field with the Navier-Stokes equa- 
tions 

d ~ * 1 
u ■ — = —u ■ L — u ■ V(p + —u ■ u) + u ■ vAu (22) 

This will simplify to 

df 

■jj- = -U- V(p + /) + u- vAu (23) 

As we know that the energy density is superharmonic on C, we deduce from the 
backwards heat equation 

|>0 (24) 
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It must be so that 

u-vAu>u- V(p + /) (25) 
The well known result is, see [5] 

■ Q < —u ■ Au + Af (26) 

Therefore we must have 

u ■ Au < Af — —to ■ io (27) 

As we know that in the set C energy density is superharmonic, that is, A/ < 
and that the enstrophy density is strictly positive in C we have 

u-Au<0 (28) 

Now the energy density evolution obtained from the Navier-Stokes equations 
must be equal to the energy density evolution deduced from the backwards heat 
equation in the set C. Therefore in the set C the the following equality must 
hold 

- u ■ V(p + f)+u- vAu = a(t)(Ap -u-Au-w-u) (29) 
Let us rearrange 

a{t)Ap = -u ■ V(p + f) + (v + a{t))u ■ Au + a(t)u ■ Co (30) 
So that we have 

g-V(p + /) , (u + <r(t)) ^ 

Ap = 1 — u-Au + lo-uj (31) 

a{t) o(t) 

We now consider the freezed system at time t = T — e and we set a(t) = a. Let 
us consider now the maximal ratio of the following in C 

max — — — s- !■ 32 

sec | u • Au | 

where u- Au ^ 0. We would not actually need the absolute value in the nomina- 
tor as it is always positive in C following from the inequalities presented earlier. 
Let the maximal ratio be 



- -. u-V(p+f) 



max ■ 



Sec I u ■ Au 



= M (33) 



As u ■ Au < is always satisfied in C, we can force the term Ap be to negative 
in C by choosing large enough value for the viscosity parameter v. This forcing 
is satisfied if 

V -±^ > M (34) 
a 

So we choose naturally 

v = oM (35) 
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Now what we have actually is that the pressure is superharmonic in C as we 
have Ap < 0. This means that the function p + f is also superharmonic in C. 
Therefore we can deduce (using Weierstrass theorem and theory of superhar- 
monic functions) that there is an interior maximum for the function p + f and 
therefore for some x S C the gradient must vanish 

V(p + /) = (36) 

So at some point in the set C at time t = T — e we must have that 

% = aMu ■ Au (37) 

Which implies that 

I < (38) 
Which is a contradiction as we had that in C at time t = T — e we have 

ft > ° (39) 

Which from we conclude that such energy cascade singularity process is impos- 
sible for the viscous, incompressible Navier-Stokes equations with such a choice 
of viscosity. 

Because the Navier-Stokes equations exhibit scale invariance, applying a renor- 
malization to the velocity field one deduces global regularity. 
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